An efficient scheme for tiie deterministic maximal entanglement of N trapped ions 
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We propose a method for generating maximally entangled states of A*' two-level trapped ions. 
The method is deterministic and independent of the number of ions in the trap. It involves a 
controUed-NOT acting simultaneously on all the ions through a dispersive interaction. We explore 
the potential application of our scheme for high precision frequency standards. 
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The entanglement of quantum states of two or more 
particles, aside from being of intrinsic interest, is of great 
practical importance in the fields of quantum cryptogra- 
phy and quantum computation [Q . One other area where 
entangled quantum states may have a significant inipact 
is that of the improvement of frequency standards P-H . 
Advances in cooling and trapping of ions have given rise 
to new techniques in high precision spectroscopy which 
may yield frequency standards with accuracies of the or- 
der of one part in 10^^-10^^ Key to the improvement 
of frequency standards beyond the shot-noise limit ||^ is 
the establishment of an entangled state of a collection 
of N two-level atoms. Initial theoretical investigations 
examined the use of squeezed spin states . We con- 
centrate here on a maximally entangled A'^-particle state, 
having the form 

|*M> = 4^ ||ei,e2,...eAr) -I- e^-^jgi, 32, •••5Jv) } , (1) 
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where \ej) and \gj) denote the excited and ground states 
of the jth particle respectively. Using the Dickc angular 
momentum states Iql this state can be written as 
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(2) 



where J = N/2. The above state is an A^-particle ver- 
sion of the Greenberger-Horne-Zeilinger state and has 
been shown to display extreme quantum entanglement 

. It may also be considered a special case of the atomic 
Schrodinger cat states 1^ . Recently, Turchette et al. have 
reported the generation of a non maximally entangled 
two-particle state using a deterministic method in an ion 
trap experiment ||l0[| . 

The maximally entangled state given in Eq. ^ may 
be used in high precision spectroscopy to measure the 
transition frequency = {E^ — Eg)/h, where and 
Eg are the respective energies of the electronic excited 
and ground states [Q . In contrast to measurements with 
uncorrelated atoms, which yield an uncertainty in the 
frequency that depends on A^~^/^, the state \^m) al- 
lows one to measure the transition frequency to an un- 
certainty of 7V~^. Huelga et al. |1^] have described how 



in principle a standard Ramsey spectroscopy scheme |T^] 
can be modified to achieve this limit. The key to this 
is a controUed-NOT (C-NOT) operation right after the 
initial and before the final Ramsey pulse, where the elec- 
tronic state of the ion which is manipulated through the 
Ramsey pulses acts as a control to flip the electronic state 
of all the other ions. 



As Huelga et al. U% have shown, in the presence of 
decoherence the standard Ramsey spectroscopy measure- 
ments on uncorrelated atoms and measurements on the 
maximally entangled states yield the same resolution. Of 
course, if decoherence is not present or if the measure- 
ments can be done in a time short compared to the deco- 
herence time, the maximally entangled states will yield 
higher resolution frequency measurements and thus it is 
of interest to find efficient mechanisms for their genera- 
tion. 

In view of recent experimental progress pC|], perhaps 
the most promising physical system for the generation of 
the type of maximally entangled state given in Eq. (||) is 
a string of laser-cooled ions in a linear rf trap. Previously, 
Cirac and ZoUer have proposed a method which se- 
quentially performs N C-NOT operations between the 
internal states of pairs of ions in such a string. This ap- 
proach requires individually addressing all the ions with a 
well focused laser beam. Bollinger et al. Q have proposed 
an alternative method that does not require interacting 
with the ions individually. It does require the use of three 
vibrational modes and the generation of linear couplings 
between pairs of those modes as electronic transitions are 
driven for all the ions simultaneously Both in the 

scheme by Cirac and Zoller [O and the method proposed 
by Bollinger et al. [Q the number of steps or laser pulses 
required to generate the maximally entangled state |^'_a/) 
is proportional to the number of ions. Only recently, 
Wineland et al. ^ have proposed for the first time a se- 
quence of operations which accomplishes this with a fixed 
number of steps. 

In this paper, wepresent a method of generating states 
of the form of Eq. (^, which is independent of the number 
of ions, and in contrast to most experiments that generate 
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entangled states, is deterministic ||l5|. While our scheme 
is related to the proposal in |^ through the sequence of 
operations which lead to the maximally entangled state, 
it is significantly different in that (i) it does not rely on a 
specific value of the Lamb-Dicke parameter, (ii) it oper- 
ates only on two electronic levels of the trapped ions and 
(iii) it points out a C-NOT operation which can operate 
on multiple ions. Our scheme requires that one of the 
ions be addressed individually for certain manipulations 
and we assume here that this is done with a well focused 
laser . We further assume that all the ions can be ad- 
dressed simultaneously with a laser beam of sufficiently 
broad waist The maximally entangled state \"^m) is 
generated through a sequence of five laser pulses, start- 
ing from an initial state |J, — J)|0), where all the ions 
are in their ground state and the collective motion has 
been cooled to the ground state . At the heart of our 
preparation scheme lies a C-NOT operation between the 
collective vibrational motion and the internal states of 
all the ions. This is generated by a type of dispersive in- 
teraction between those degrees of freedom. Previously, 
this kind of coupling has been discussed in connection 
with the generation of the vibrational Schrodinger cat 
states for a single trapped ion and as a degenerate 
Raman-coupled model in the context of cavity quantum 
electrodynamics pOj . 

Before going through the preparation scheme step by 
step, we consider the four types of pulses involved. We 
require pulses at the two frequencies wq and loq — i^x, 
where is the frequency of the ions' collective harmonic 
motion along the trap axis which we take to be the x 
axis We assume that all pulses are performed with 
laser beams that are derived from the same source, so 
that at time t — 0, all electric fields that excite the 
system during the preparation scheme are in phase. To 
derive the transformations caused by the various pulses 
we generically assume the respective laser beams to be 
turned on from time t = to, to t = to + tp, and we give 
the corresponding unitary transformations U{to,tp), in 
a rotating frame where = exp[iHfjt/h]\'if{t)), 

Hr — ?it^o X]i3=i and \^{t)) is the state of the 

system in the Schrodinger picture. In our analysis of 
the dynamics generated by the laser excitation we ex- 
plicitely consider only the collective motion along the 
trap axis, characterized in terms of the vibrational en- 
ergy eigenstates \n). The motional state perpendicular 
to the trap axis remains in the ground state throughout 
the preparation scheme since we do not excite sidebands 
of the motion along those directions and in a linear trap 
Vx <C jAi , I'z ■ We further assume both the Lamb-Dicke 
limit |2^] and the low excitation regime |^^. The first 
condition allows us to expand the Hamiltonian describ- 
ing the interaction of the ions with the laser light in terms 
of the Lamb-Dicke parameter rjx = fc^Azo, where kx is 
the projection of the laser wave vector onto the trap axis 
and Acco = {Ti/2Mi'xY^'^ , is the width of the motional 
ground state for a single ion along that axis, M being 
the mass of a single ion. In the low excitation regime 
we retain only resonant transitions in the analysis of the 
excitation. 



We first consider the action of the pulses involving only 
one of the ions. The first of these is a resonant 7r/2-pulse 
where a laser beam of frequency loq, propagating perpen- 
dicular to the trap axis excites, say, the A^th ion. This 
generates the transformation 

oo 
n— 

® {lew) (bat I -I- \gN){gN\ 

- l5Jv)(e7v| + \eN){9N\} , (3) 

where the pulse duration — i^ 1'^^, and U, denotes the 
Rabi frequency for that pulse. The phase factor is due 
to the free evolution of the vibrational degree of freedom 
during the pulse. 

The second type of pulse required is a 7r-pulse gen- 
erated by a laser of frequency ujq — Vx, and where the 
wave vector has a component kx along the trap axis. 
This generates a \gj^)\n -\- 1) ^ \eN)\n) transition of 
the Jaynes-Cummings type, where the Rabi frequency 
rij^ oc rix\/n+ l/\fN . For a given \gM) \n+\) <^ |eAr) |n) 
transition a pulse of duration t'')^ ^ = 7r/f2jp, causes the 
transformation 

+ i e""-^ V°-^'''cA \n){n + l\® \eN) {gN \ , (4) 

where the phase factors are due to the free evolution of 
the vibrational degree of freedom during the pulse and 
the fact that we have assumed all electric fields to be in 
phase at time t = 0. We note that the state Isat) |0) is not 
coupled by the Jaynes-Cummings pulse. 

The third type of pulse which we shall refer to as a 
dispersive 7r-pulse is generated by two laser beams of fre- 
quency loq . More specifically, the first beam is propa- 
gating perpendicular to the trap axis and the second one 
has a wave vector component kx along that axis. While 
not exciting any vibrational sidebands the pulse exploits 
the dependence of the generated Rabi oscillations be- 
tween the states |5Ar)|'^), and |eAr)|n), on the motional 
excitation number n, which arises from the spatial vari- 
ation of the electric field along the x axis As shown 
in [|l8|, if the two laser beams have a relative phase dif- 
ference of TT their amplitudes can be chosen such that 
the (spatially) constant terms of the electric fields as- 
sociated with the two laser beams cancel each other and 
the Rabi frequency fl^J^ °^ Vx''^/^^ to leading order in the 
Lamb-Dicke parameter . In deriving the pulse trans- 
formation we have assumed here the second laser beam 
to be phase shifted by tt, with respect to the first. For 
a given \gN)\n) 4=> |eAr)|n) transition a pulse of duration 

^'dil TT ~ ^/^'dil' then generates the transformation 

{|ejv)(.gAr| - \gN){eN\} , (5) 
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which is effectively a C-NOT operation between the col- 
lective motion being in either of the states |0) and \n), 
and the electronic state of the A^th ion. For a motional 
state |0), the Rabi frequency ri[j°] = 0, and the elec- 
tronic state remains unaffected. For the motional state 
\n), the above transformation flips the electronic excited 
and ground state, apart from a phase factor. We note as 
an important feature of the C-NOT operation proposed 
here that it remains valid even beyond the Lamb-Dicke 
limit since it's essential feature of not affecting the mo- 
tional ground state |0) does not depend on the exact value 
of the Lamb-Dicke parameter. 

The fourth type of pulse required in our preparation 
scheme is a dispersive 7r-pulse as described above, but 
acting on all the ions simultaneously. The set-up of the 
two laser beams that generate this pulse is identical to 
the single ion case described above but where the beam 
waists are assumed sufflciently broad to excite all the ions 
with the same strength . We further assume that the 
spatial variation of the exciting electric fields along the 
trap axis is small compared to the separation of the ions 
in the trap, i.e. fcajA^, <C 1, where the separation A^;, is 
typically several /xm in current trapped ion experiments 
|^,|6|. This may be realized by having the second laser 
beam propagating almost perpendicular to the trap axis 
or through a Raman excitation which allows to control 
the effective wave vector . Alternatively, one can con- 
sider using two beams whose wave vectors are sensitive 
to the two (independent) radial directions and employing 
the ions' collective motion perpendicular to the trap axis 
p6t . In either situation the phase of the exciting electric 
fields is equal for all ions |^ , and the dynamics generated 
by the pulse is given by the product of the single-ion time 
evolution operator given in Eq. (^). The dispersive inter- 
action then couples the states \ J, —J)\n) , and \ J,J)\n), 
and for a given transition the transformation 



in) N 
Ndis,7T y 



\n){n\ 



N 



(6) 



is generated through a pulse duration t^jj^ ^ — T^/^^Ndts^ 



where 



in) 
Ndis 



In/N, denotes the Rabi frequency for 



that transition 

We now go through our preparation scheme step by 
step, starting from the initial state |4'ii;(0)) = | J, — J)|0). 
Since the scheme involves acting on one of the ions 
separately, we shall, when necessary write |J, — J) = 
\J' ,—J')\gN), where J' = {N — l)/2. First, a resonant 
7r/2-pulse is applied to the iVth ion to produce the state 



|vI/«(ii)) = ^|J', 



-^'>{M + |ejv)}|0) 



(7) 



at the time ti = Next, a Jaynes-Cummings 7r-pulse 
transfers the superposition of the electronic state of the 
A^th ion into the collective motion along the trap axis. 
From Eq. (^) with n = 0, the resulting state at time 

t2 = h 



t 



(0) 

JC,7ri 



\^R{t2))^^\j,-j){n 



(8) 



Then, the superposition of vibrational states is trans- 
ferred into the electronic degrees of freedom of all the 
ions simultaneously by applying a dispersive 7r-pulse to 
all the ions. As seen from the pulse transformation 



^wis(^2,i';v^is^^), given in Eq. (|^) this effectuates a C- 
NOT operation between the collective vibrational state 
and the internal states of all the ions simultaneously, in- 
dependent of the number of ions. The state after the 
pulse is 



1 



|vI'^(t3))-^{|J,-J)|0)+ze- 



\J,J)\l)]. (9) 



where ts = <2 + ^ivdis tt ■ have now generated the re- 
quired superposition between the Dicke states |J, — J) 
and I J, J) . The remaining two pulses serve to disentangle 
the vibrational and electronic degrees of freedom. The 
first of those is a dispersive 7r-pulse acting on the A^th 
ion. From Eq. (^) with n = 1, the state resulting from 

this pulse at time ti~t^ + t^^^^ ^, is given by 



\^,j,it,))^—{\j',-j')\g^)\0) 

-»e— ^*^|J',J')|.9^)|1)} 



(10) 



Finally a second Jaynes-Cummings 7r-pulse, identical to 
the one that led us from Eq. (^ to Eq. (||) , realizes the 
maximally entangled state 



\^R{t,)) = —{\J,-J) + \J,J)}\0). 



at time t^ = t^+t^j^ ^, leaving the collective motion along 
the trap axis in its ground state. In the Schrodingcr pic- 
ture |*(t5)) = |*A/), where cj) = iVwo^s, in Eq. (|). The 
phase (f) in the maximally entangled state can be con- 
trolled by changing the phase of the initial 7r/2-pulse 
with respect to the other electric fields. Note that the 
maximally entangled state is produced deterministically 
by the procedure described here. Moreover, all pulses 
considered here drive the same |e) <^ \g) transition. This 
is important for the experimental realization of our pro- 
posal since single transitions can be made independent 
of magnetic field fluctuations to first order . 

As we have said earlier, the maximally entangled state 
generated in Eq. (|l^) may be used in high precision spec- 
troscopy |^,Il[. In Ref. BoUinger et al. describe a 
Ramsey technique where once the maximally entangled 
state has been established two Ramsey pulses are applied 
to all ions simultaneously and the expectation value of 

the product operator Hil:! — \9i){9i\} i serves to 

extract the transition frequency loq. This is measured by 
determining the number of ions in the excited or ground 
states. In order to not degrade the signal-to-noise ra- 
tio, the uncertainty in this measurement must be <C 1 



(11) 
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atom which requires that the number of ions in the trap 
be small. In contrast, the Ramsey technique described 
by Huelga et al. relies on population measurements 
on a single ion and the state generation is an integral 
part of the modified Ramsey scheme. In this situation 
the role of the initial 7r/2-pulse at frequency luq is taken 
over by a 7r/2-Ramsey pulse at frequency to, which is 
detuned by a small amount A — luq — to, from the res- 
onance frequency loq which one aims to determine. We 
therefore assume that the pulses which prepare the max- 
imally entangled state are generated by laser beams at 
the frequencies uj and uj — v^. Thus the frequency wq is 
replaced by to in the pulses, and we describe the time 
evolution in a frame rotating with the frequency lo of 
the Ramsey pulses. We denote the state of the sys- 
tem in this frame by and under the assumption 

|A| <C il, rij*^, ri^j"], ri^jjj,, the pulse transformations in 
that frame are the same as given in Eqs.(||)-(^. Start- 
ing from the state |^_r'(0)) = \J,—J)\0) the pulse se- 
quence described above leads to the state |4'fl'(t5)) = 
{I J, - J) + I J, J)} \0)/V2. We then assume that the sys- 
tem is let free to evolve for a time T, resulting in 
the state |*fl'(t6)) = {\J,-J)+e-'^^^\J,J)}\0)/V2, 
where — tr^ + T. Then the pulse sequence which gener- 
ates the maximally entangled state is applied again but 
in reverse order. This results in the final state 

+ (l-(-ire-^^^)|eA,)}|0), (12) 

at the time = T + 2^5. From this state the reso- 
nance frequency ujq can be determined with uncertainty 
6ujQ cx 1 /iV, , by measuring the internal state of the 
mh ion, and where P = {1- (-l)^cos [iVAT]}/2, gives 
the probability of finding the A^th ion in its excited state. 

To address the susceptibility of our scheme to imperfec- 
tions in the state generation procedure we emphasize that 
the experiment reported in |]l0|] has demonstrated the 
generation of almost maximally entangled states through 
unitary manipulations. The observed fidelities of approx- 
imately 0.7 indicate that deviations from unitary evolu- 
tion do not destroy the sought-after state and affirm the 
validity of the unitary analysis presented here. 

In summary, we have proposed an efficient method for 
generating maximally entangled internal states of a sys- 
tem of N trapped ions. The method has the further ad- 
vantage of being deterministic. Finally, we have shown 
how a such a state and the generation scheme described 
here may be used for high precision Ramsey spectroscopy. 

ACKNOWLEDGEMENTS 

This work was supported in part by the UK Engineer- 
ing and Physical Sciences Research Council and the Eu- 
ropean Community. C.C.G. acknowledges the support of 
a PSC-CUNY grant and a Schuster Award from Lehman 
College. J.S. is supported by the German Academic Ex- 
change Service (DAAD-Doktorandenstipendium aus Mit- 
teln des dritten Hochschulsonderprogramms). We wish 



to thank Drs. S.F. Huelga, M.B. Plenio, P.L. Knight and 
D.J. Wineland for helpful discussions. 



[1] A. Barenco, Contemp. Phys. 37, 375 (1996); R.J. Hughes 
et al, ibid, 36, 149 (1995); M.B. Plenio and V. Vedral, 
ibid, to be published (1998). 

[2] D.J. Wineland et al, Phys. Rev. A 46, R6797 (1992); 
W.M. Itano et al, ibid, 47, 3554 (1993); D.J. Wineland 
et al., ibid, 50, 67 (1994); D.J. Wineland et al., e-print 
quant-ph/980902^ . 

[3] D.J. Wineland et al. Journal of Research of NIST, 103, 
259 (1998). 

[4] J.J. Bollinger et al, Phys. Rev. A, 54, R4649 (1996). 

[5] D.J. Wineland et al, IEEE Transactions on Ultrason- 
ics, Ferroelectrics and Frequency Control, 37, 515 (1990). 
D.J. Berkeland et al, Phys. Rev. Lett. 80, 2089 (1998). 

[6] R.H. Dicke, Phys. Rev. 93, 99 (1954). 

[7] D.M. Greenberger, M.A. Horne, and A. Zeilinger, in 
Bell's Theorem, Quantum Theory, and Conceptions of 
the Universe, edited by M. Kafatos (Kluwer, Dordrecht, 
1989); D.M. Greenberger, M.A. Horne, A. Shimony, and 
A. Zeilinger, Am. J. Phys. 58, 1131 (1990). 

[8] N.D. Mermin, Phys. Rev. Lett. 65, 1838 (1990). 

[9] C.C. Gerry and R. Grobe, Phys. Rev. A 56, 2390 (1997); 
ibid, 57, 2247 (1998); G.S. Agarwal, R.R. Puri, and 
R.P Singh, ibid, 56, 2249 (1997); M.G. Benedict, et al. 
Acta Phys. Slovaka 47, 259 (1997). 
[10] Q.A. Turchette et al, Phys. Rev. Lett. 81, 3631 (1998). 
[11] S.F. Huelga et al, Phys. Rev. Lett. 79, 3865 (1997). 
[12] N.F. Ramsey, Molecular Beams, (Oxford University, Lon- 
don, 1963). 

[13] J.I. Cirac and P. ZoUer, Phys. Rev. Lett. 74, 4091 (1995). 
[14] J. Steinbach, J. Twamley, and P.L. Knight, Phys. Rev. 

A 56^815 (1997). 
[15] See |lC|| and references therein. 

[16] R.J. Hughes et al, Fortschritte der Physik 46, 329 
(1998). For the two particle case, Turchette et al. |ic| ] 
have experimentally demonstrated an alternative tech- 
nique to differentially address single ions by controlling 
the ion micromotion. 

[17] B.E. King et al, Phys. Rev. Lett. 81, 1525 (1998). 

[18] C.C. Gerry, Phys. Rev. A 55, 2478 (1997). 

[19] P.L. Knight, Phys. Scr. T 12, 51 (1986); S.J.D. Phoenix 
and P.L. Knight, J. Opt. Soc. Am B 1, 116 (1990). 

[20] We neglect the ions' micromotion whose main effect is the 
alteration of transition rates between quantum levels |2^ 
which can be accounted for by experimental calibration 

i- 

[21] P.J. Bardroff et al, Phys. Rev. Lett. 77, 2198 (1996). 

[22] W.E. Lamb, Phys. Rev. 51, 187 (1937); R.H. Dicke, 
Phys. Rev. 89, 472 (1953). 

[23] J.I. Cirac et al. Adv. At. Mol. Opt. Phys. 37, 237 (1996). 

[24] The Rabi frequency here is reduced through its depen- 
dence on rj^/N, which may be compensated for by high 
enough laser power. Second sideband excitation has been 
demonstrated for a single ion in Ref. [ p5[ . 

[25] D.M. Meekhof et al, Phys. Rev. Lett. 76,1796 (1996). 

[26] C.R. Monroe and D.J. Wineland, private communication. 



4 



